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Critical exponents of the Kawasaki dynamics in the Ising chain are re-examined numerically 
through the spectrum gap of evolution operators constructed both in spin and domain wall repre- 
sentations. At low temperature regimes the latter provides a rapid finite-size convergence to these 
exponents, which tend to z ~ 3.11 for instant quenches under ferromagnetic couplings, while ap- 
proaching to z ~ 2 in the antiferro case. The spin representation complements the evaluation of 
dynamic exponents at higher temperature scales, where the kinetics still remains slow. 
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I. INTRODUCTION 

Kinetic Ising models have long been helpful in understanding general questions of nonequilibrium statistical physics, 
as well as for elucidating particular issues of experimental relevance There is by now a vast body of research 
studying these models under Glauber and Kawasaki dynamics @, Q , the so called models A and B respectively in 
the terminology of Hohenberg and Halperin The first type of dynamic considers single spin flip processes for 
describing relaxation towards equilibrium in a variety of magnetic materials, while the second case consists of spin 
exchanges specially aimed to study stochastic processes under conservation of total magnetization. This constraint 
makes relaxation slower, and is instrumental in studying phase separation, domain growth, and freezing observed 
after rapid cooling in several systems, such as binary fluids and alloys All these phenomena are characterized by 
a long relaxation time r which near critical points diverges as r oc £ z , where £ denotes the spatial correlation length 
(proportional to the typical system size), and z is the dynamic exponent of the universality class to which the dynamic 
belongs. 

When it comes to one dimension (Id), these dynamics are still largely amenable to experimental probe 5] apart from 
being interesting in their own right Q , and as is well known, in the Glauber case allow for an exact treatment 0, H| ■ 
In contrast to that latter aspect, the equations of motion for the Kawasaki dynamics conform an infinite hierarchy 
which to our knowledge do not admit analytic solution, so the problem can only be approached from approximations 
PJ. Even though there is no finite critical temperature through which to quench, using random walk arguments to 
derive an approximate kinetic equation for the kink density under ferromagnetic (F) interactions [7|, interestingly, 
the dynamic exponent has been put forward to be the same as observed in higher dimensions, namely z = 3, after a 
deep quench to a low temperature (LT) . Although for d > 2 such exponent was also accounted for by renormalization 
group arguments [l|, Q as well as by surface dynamical considerations @, other studies in Id related systems 
suggest instead a value of z = 5 ll], in agreement with linear response schemes 

In this work we re examine this exponent by constructing and diagonalizing numerically the kink evolution operator 
associated to the master equation of the Kawasaki dynamics in finite chains. Following the physical picture given in 
0, at late evolution stages and LT regimes such operator would essentially describe a highly diluted system as the 
kink or domain wall density is basically a measure of the average inverse domain size. Consequently, in evaluating 
dynamic exponents from finite samples it might be expected that size effects will pose no severe limitations for the 
kink representation. In fact, as we shall see, the kink density will get fairly small even for the low lying levels of 
the evolution operator, that is what ultimately matters at large times. On the other hand there is also the question 
about estimating dynamic exponents at finite temperatures for which the Kawasaki kinetics still remains critical, as 
opposed to the exponential decay of the Glauber dynamics in Id. Although in this situation dynamic exponents are 
no longer related to domain growth - these are cut off by a finite correlation length - they still provide the fundamental 
relation between the typical system size and its relaxation time. We shall also address this issue by diagonalizing the 
evolution operator of the original spins. Insofar temperature is not too low, these latter become weakly correlated 
and finite-size effects will not be paramount. A posteriori, our results will lend further support to this view. 

Therefore these dual representations - kinks and spins - provide a means for probing the robustness of dynamic 
exponents throughout different regimes. As a result, it will turn out that for F couplings z varies continuously from 
~ 3.11 to ~ 2 at low and high temperature (HT) scales respectively. Such nonuniversality should come as no surprise, 
since in the HT limit the Kawasaki dynamics simply reduces to a diffusive disordering kinetics. However under 
antiferromagnetic (AF) exchanges these exponents are rather robust and, in line with previous findings jlOj , their 
values will remain diffusive even at LT regimes, just as the z of the Glauber universality class. We shall return to this 
question later on within a context of antikink operators. 

The layout of this work is organized as follows. In Sec. II we recast the Id- Kawasaki dynamics in terms of a quantum 
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spin analogy that readily lends itself to evaluate spectrum gaps of evolution operators, i.e. dynamic exponents, both 
in the spin and domain wall representations. Due to detailed balance, either of these descriptions can be brought to a 
symmetric representation by means of simple non-unitary spin rotations. This simplifies considerably the subsequent 
numerical analysis of Sec. Ill in which spectrum gaps are obtained via standard recursive techniques (l3| in various 
situations. For now let us remark that already modest chain lengths are able to yield clear finite-size trends both at 
LT and HT scales. Finally, Sec. IV contains a summarizing discussion along with some remarks on extensions of this 
work. 



II. SPIN AND KINK REPRESENTATIONS 



As is well known, the Ising model has no intrinsic dynamics because all spin operators involved in its Hamiltonian 
commute with one another. Therefore, an ad hoc dynamics must be prescribed by coupling the system to a heath 
bath at temperature T so as to induce energy changes in the model. This is described by a master equation [hH ] for 
which we introduce briefly here some preliminary considerations. 

Basically, the dynamics is associated to a gain-loss equation, constructed generically as 

8 t P(s, t) = [ W(s' -> a) P(s', t) - W(s -> a') P(a, t) } , (1) 

s f 

which governs the evolution of the probability P(s,t) that the system will be at state \s) at time t. The elementary 
change steps are embodied in the transition probability rates W(s — > s') per unit time at which configuration \s) 
evolves to \s') . For our purposes, it is convenient to think of this equation as a Schrodinger evolution in an imaginary 
time, namely d t \P{t) ) = —H\P{t) ), under a pseudo Hamiltonian or evolution operator H = H d + H n d whose diagonal 
and non-diagonal matrix elements are given by 

(a\H d \a) = J2 W(a^a'), (2) 

( s' | H nd \s) = -W{s^s'). (3) 

Formally, this enables one to derive all subsequent probability distributions \P{t) ) = J2 S P( s i I s ) from the action of 
H on a given initial distribution, i.e. \P(t) ) — e~ H *|P(0) } [1J]. In particular, the relaxation time of any observable 
is singled out by the eigenvalue \\ corresponding to the first excitation mode of H, i.e. r _1 = ReA x > 0, whereas 
the steady state merely corresponds to an eigenvalue Ao = [l4j. If the steady configuration should actually 
coincide with the Boltzmann equilibrium distribution, the above rates must be constrained by detailed balance, that 
is W(s — > s')e~^ E ^ = W(s' —> s)e~P E ( s Vs, s', where E stands for the respective energy configurations of the 
system in question, and f3 = 1/ksT. 

Detailed balance itself can not determine entirely the form of such rates, thus for the specific case of the Kawasaki 
dynamics hereafter considered we take up the common choice [l|, |3| 
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W(s -»■ a') = - 1 - tanh £ AE S ^ , (4) 



where \s), \s') are states differing at most in a pair of nearest neighbor (NN) spins exchanged at some location 
i = 1, 2, ... L, and AE S , S / = E(s') — E(s) is the change of the Ising energy E = — J J^i s i s i+i> ( s i — ±1) either with 
F (J > 0), or AF (J < 0) interactions. Therefore depending on the spin states at locations i — 1 and i + 2, the rate 
at which spins (s^, s i+1 ) = (s, — s) exchange their states results in 

«.-»«(-) = i^(— -.), 

where K = /3J, a K = tanh2A', and the signs ± denote forward and backward hoppings as depicted schematically in 
Fig. 1. If we think of the spin configurations \s) as being already diagonal in the z-direction, say, then by promoting 
Si variables to Pauli matrices of , the operational analog of Eq. (j3]) will read 

H nd = - E [ R i-hi+2 (K) + RU,+2{-K) o-+<t~ +1 ] , (6) 

i 

where <7 + ,a~ are the usual spin-i raising and lowering operators, while Rf_ 1 i+2 (±A") = 5 ± -^ L (cf +2 — o-f_ x ), 
is simply the operational counterpart of Eq. ([5]). As for the diagonal elements of © needed for conservation of 
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Ri-l,i+2( K ) 

Ri-u+2(- K ) 



FIG. 1: Transition probability rates for neighboring spin exchanges, as defined in Eq. 



probability, notice that these basically count the number of hoppings in which a given configuration can evolve to 
different ones by exchanging NN spins at a time. This can be properly tracked down in terms of number operators 
h = tj + o~ probing and weighting all NN spins, namely 

ff«i = E[^-M+aW^( 1 -^i) + ^-i,^2(--«0*W+i(l-^)] • (7) 

i 

Although the correlated exchange terms of Eq. ^ leave us with a non-symmetric evolution operator, in preparation 
for the numerical analysis of Sec. Ill we can make some progress by exploiting detailed balance. This latter warrants 



the existence of representations in which H is symmetric and thereby diagonalizablc 14]. For our purposes, it suffices 
to consider the diagonal non-unitary similarity transformation 

S = e^E^W+i j ( g ) 

under which the hopping terms of H n( i transform as 

af +l af -> exp [ (af +2 - af_j) ] af +l af , (9) 

while leaving invariant all the above mentioned diagonal operators. Hence, after straightforward manipulations it 
can be readily verified that H SJ3 j n = SHS^ 1 actually produces a self adjoint spin representation which for periodic 
boundary conditions (PBC) is found to be 

Hapin - -i(l+7 K )^ (l+tanh 2 ifaf_ 1 af +2 )( CT ? CT ? +1 +aX+i) 

i 

+ 4 12 [ 1 + a * «+2 - ( X + "kX^+i ] > ( 10 ) 

i 

where j K = sech2i4T. In passing, it is worth pointing out that in the HT limit this operator coalesces into a 
Heisenberg ferromagnet, i.e. H sp in — > —\ Ylt (&i ' ^t+i — 1)> with a spectrum gap behaving like cx 1/L 2 , as it should 
for a disordering diffusive kinetics. On the other hand, in such regime the equilibrium correlation length becomes 
of the order of the lattice spacing, thus the numerical diagonalization of H sp in in finite chains can be a reasonable 
starting point to evaluate dynamic exponents as temperature is quenched slightly (though instantly to a constant 
value, because transition rates have been taken time independent throughout). 

When it comes to LT regimes however, Eq. (|10[) is not of much practical use, since the system approaches a state 
consisting of long magnetic domains whose typical sizes are of the order of the correlation length £ ~ e 2 '^' (nonetheless, 
for J > see the upper bounds to z and their extrapolations provided in Sec. Ill A). In this more interesting scenario 
the dynamics is strongly affected by whether the coupling exchanges are F or AF, thus it is convenient at this point 
to separate the discussion accordingly. 



A. Kink representation: J > 



As was mentioned in Ref. [7] , the dynamics of this situation is basically mediated by single spins detaching from a 
domain wall (a kink), and then performing a random walk. Some of these spins will return to their original domains, 
while others will reach the next ones. On average, as a result of many of these processes the domains themselves 
will perform a random walk until they meet and merge into larger domains. Besides, single meandering spins may 
eventually nucleate and trigger the growth of a new domain. The key issue in this description is that for large times 
very few kinks will survive at LT scales. In particular, in equilibrium these are totally uncorrelated and with an 
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average density ~ l/£, whereas that of the spins mediating the whole process is of order l/£ 2 as they are made 
up of two consecutive kinks (see Fig. [5^ below). Thus, for large times it makes sense to attempt a nonequilibrium 
description in terms of kinks rather than spins. Presumably (as confirmed later on in Sec. Ill), the first excitation 
mode of the kink evolution operator will also correspond to a diluted state, so size effects in a numerical estimation 
of z will not be as severe as in the spin description. 

1/2 
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FIG. 2: Transition rates for (a) diffusion of double kinks, and (b) assisted deposition-evaporation of kinks. These latter are 
denoted by vertical lines separating domains of opposite spin orientations. 

To recover the equilibrium behavior of the Ising model from now on we consider the case of zero magnetization, 
so this dual representation corresponds to the two-to-one mapping outlined schematically in Fig.[2j If we imagine 
these kinks as hard core particles A, then the Kawasaki dynamics involve two basic processes (i) dimer diffusion 
A + A + $^$ + A + A, and (ii) 'assisted' deposition-evaporation ft + A + fl^A + A + A, as schematized respectively 
in Figs. [5^ and[5J). The first situation represents the meandering spins referred to above and involves no energy changes 
(rates 1/2). Notice that these dimers (double kinks) do not preserve their identity as they may eventually dissociate 
when contacting a domain wall (single kink). The second process corresponds to a spin detachment (attachment) 
from a domain wall, and its rates e K (e_ K ) parallel those of Eq. ([5]) when S{_i = — Si+2, thus 



e±K = - ( I =p tanh 2K ) . (11) 
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Therefore, following the reasoning steps discussed for the spin representation and reinterpreting kinks as ^-spinors 
in a dual chain, clearly the non-diagonal terms H® d + H^ d of the evolution operator associated to the (a) and (b) 
processes of Fig. [2] can be constructed as 

H a nd = ~ J> ( 4+x^x + h-c.) , (12) 

i 

H nd = ( C K a i-l a i+l + € -K a i-l U i+\) > ( 13 ) 

i 

with number operators acting here as projectors that rule out vacancy mediated processes. To recast the non- 
diagonal operator into a symmetric representation we recur once more to detailed balance and rotate all 's around 
the z direction using a common pure imaginary angle ip = iK . This is carried out by the non-unitary similarity 
transformation 

S = ^ °* . (14) 

for which it is straightforward to show that 

4-1^+1 ->• e ±2 V-i<i , (15) 

thus producing the symmetrization of H^ d while keeping all terms of H^ d (already hermitian) unaltered. 

As before, preservation of probability is taken into account by a diagonal operator H d + H d balancing both of the 
above (mutually exclusive) events, in turn probed respectively by 

H d = 7j [nj-i fjj (1 - h i+ x) + (1 - hj-x) hj n i+ i ] , (16) 

i 

H\ = 22 [ e K (1 - fk-x) hi (1 - h i+1 ) + e_ K fk-i hi h i+1 ] . (17) 
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Since the imaginary rotation (fl4l) is also diagonal, it has no effect on these latter operators. Thus, collecting all terms 
and using PBC throughout, after some algebraic steps we finally obtain a self adjoint representation for the kink 
evolution operator, namely 

n kink = -^(l + o?) [(l + 7x)of-iof+i + (l-7ir)oti^i] 

i 

+ iE( 1+a --<<+i) + ^ e ---E< . K>0 - ( 18 ) 

i i 

As can be readily verified from this equation, in the limit T — > + the action of Hfa'nfc on any configuration having 
non NN kinks vanishes like 0(e~ 4K ), thus yielding a metastable state. In the spin representation this corresponds to 
configurations with domains lengths larger than the lattice spacing, which amounts to hindering spins to detach and 
so diffuse through. In turn, from Eq. (|10p it can be checked that this situation also yields metastable states of H sp i n 
so long as J > 0. Clearly, the number of these configurations grows exponentially with the system size, and most of 
them take over the asymptotic dynamics with long lifetimes oc £ 2 . This marks an important difference with respect 
to the AF dynamics to be introduced briefly in Sec. II B. 

As an aside, we finally mention that by construction Hkink n °t only preserves the parity of kinks e i7r ^-'i " J (being 
even for PBC), but in turn satisfies 
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0, (19) 



which simply expresses the conservation of the total spin magnetization in the original system. In practice, for the 
numerical evaluation of gaps we will just build up the adequate kink states from the corresponding spin ones. 



B. Antikinks: J < 



While most of the above ideas and procedures applies to the AF dynamics as well, ultimately the relaxation to 
equilibrium becomes faster than in the F case. By changing from a description based on kinks to AF domain walls, i.e. 
to antikinks, we can easily construct an evolution operator that for large times and LT regimes essentially describes 
a diluted antikink system. As before, if we think of antikinks as hard core particles A (former vacancies under F 
couplings), the basic processes now involve (i) next NN hopping, i.e. A + + 0^0 + + ^4, and (ii) deposition - 
evaporation 'assisted' by vacancies, that is + + 0^^ + + A. 

1/2 
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FIG. 3: Transition rates for (a) next nearest neighbor antikink hoppings mediated by vacancies, and (b) assisted deposition- 
evaporation of antikinks. As in Fig. 2, these latter are represented by vertical lines but separating antiferro domains. 

With the aid of these schematic events, illustrated respectively in Figs. [3^ and|3jD, and after carrying out the rotation 
referred to in Eq. (|14[) but with an argument (p = —iK, the symmetric representation of the antikink operator can be 
finally cast as 

R antl = ^(l-^f) [(1 + 7k )^-i< + i + (1-7 k )^_ 1 <i] 

i 

+ iY,( 1+a « a i a i+i)-l'-«X< a * ' K<0 - (20) 

i i 
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Of course, by exchanging the role of particles and vacancies, i.e. using the canonical transformation — > , one 
recovers the form of Eq. ([T8"]) , but this likelihood is only superficial. Note that for K — > — oo the effect of the uniform 
field e_ K vanishes like 0(e 4K ), so the only jammed state on which the action of the diagonal terms cancels out at this 
order is the vacuum of antikinks. Thus, in the limit T — > + the dynamics is no longer dominated by metastable states 
(instead proliferating exponentially for J > 0). In this situation AF domains can grow unhindered, as single antikinks 
now detach and diffuse freely (Fig. [3]), thus allowing the dynamics to run smoothly towards equilibrium. As a result 
the dynamic exponent, whose numerical analysis we next turn to consider, becomes smaller (actually diffusive) than 
for F exchanges. 



III. NUMERICAL RESULTS 



We now investigate numerically the spectral gaps of the Kawasaki operators given in Eqs. (ITUl) . (p"5|) . and (f2U|) 
along with their dynamic exponents in different situations. As a consistency check first we verified that the (rotated) 
Boltzmann distribution corresponds in fact to the steady state of those operators. This also served to start up 
the Lanczos recursion with a random state but chosen orthogonal to that equilibrium configuration. Thereafter we 
obtained the first excited eigenmodes of our symmetric H's using periodic chains of up to 24 sites, the main limitation 
for this being the exponential growth of the space dimensionalities. Another restrictive issue is that below ksT/\J\ ~ 
0.1 - 0.2 the Lanczos convergence slows down progressively because in most situations the spacing of low lying levels 
turns out to decay as e~ 4 '^'. Thus, in what follows we content ourselves with giving results above that region where 
nonetheless clear saturated tendencies can be already obtained. 



A. J > 



Turning to the F dynamics, in the insets of Figs.|4ji and [4]d we exhibit the finite-size behavior of spectral gaps 
for the case of spins and kinks respectively. They are all consistent with a gap vanishing like 1/L Z , although with 
a nonuniversal temperature dependent dynamic exponent z(T). For the purpose of observing in more detail the 
trend of size effects on these exponents under different temperature regimes, we considered a sequence of effective 
approximants defined as 

_ ln[X 1 (L-2)/X 1 (L)] 
L ~ In [L/(L- 2)] ' [ } 

(L even), which simply provide successive measures of the gap closing in either of the above representations. These 
are shown in the main panels of Figs.Hk. and 0b. As expected, in the HT region H sp j„ is able to yield fairly convergent 
estimations of z (~ 1.99), so recapturing the plain diffusive limit above T/J ~ 3 (henceforth the Boltzmann constant 
ks is set equal to 1). Besides, in this representation the above approximants turn out to constitute upper bounds of this 
exponent for almost all temperatures studied, though exhibiting different spreadings. We attempted to extrapolate 
these bounds using several procedures, but owing to the variable spreading this resulted in a rather noisy limiting 
curve. Nevertheless, in the region that most interests us the spreading becomes more stable, thus upon using a van 
den Broeck - Schwartz extrapolation scheme [l5[ we found that z ~ 3.1(0), in reasonable agreement with analysis and 
Monte Carlo simulations under instantaneous quenches Q- Moreover, this result is also consistent with estimations of 
z arising from Hkink approximants around the same region (T/ J < 0.4) where, as discussed above, the kink description 
in finite sizes is most reliable. In fact, no extrapolations are needed here since already the values of Z20 ^ 222 z 24 
follow very closely one another, collapsing near a value of z ~ 3.11. This is the main result of this section. 

As an added bonus, note that in the LT regime the kink approximants come out conforming a sequence of lower 
bounds which conveniently complement the upper ones provided by the spin representation. This is illustrated in 
Fig. 2]: where for T/J < 0.4 both kink and spin approximants should actually enclose the value of z in the limit 
L — y 00 (presumably much closer to the kink border). Using the results of our maximum reachable sizes, this suggests 
that in the ordering regime the dynamic exponent is bounded as 3.11 < z < 3.18. For 0.5 <■ T/J < 2 the trend of 
finite sizes of the kink approximants now reverses and, alike the spin Zl's, turn out to yield upper bounds of z. The 
kink approximants converge slightly faster in this region (see inset of Fig. 0b), until for T/J > 2 the trend of bounds 
reverses once more. Hence, at HT scales the actual value of z is enclosed again by our kink and spin approximants 
though this time these latter take the lead, as it should. 

As mentioned earlier, the LT convergence of kink approximants was somehow foreseeable on the basis of equilibrium 
considerations but it is not yet clear whether these apply to nonequilibrium as well. In Fig. [5] however we show that 
this also the case, at least for the sizes at hand. There, we display the kink density p x = ~ X^(V'il^l^i) for the first 
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FIG. 4: (Color online) Effective dynamic exponents of Eq. (|21[) for J > using 12 < L < 24 with L even, computed for (a) 
spins, and (b) kinks. Sizes increase downwards almost throughout (a), and upwards in (b) at both leftmost and rightmost 
temperature regions. The respective insets display the gap behavior with the system size for T/J = 10, 1, 1/2 and 1/4 (top to 
bottom). To compare slopes, data have been shifted upwards with respect to T/J = 10. Straight lines are fitted with slopes 
— Z2a{T) as calculated from the main panels. In (c) we compare these exponents for kinks and spins (dashed and solid lines). 
At low as well as at high temperatures zt and z 3 yield respectively lower and upper bounds for the thermodynamic limit of z. 
The inset of (c) provides a measure of convergence for these quantities by depicting Alz = zl-2 — zl for our maximum sizes. 
Horizontal short doted line in (a) stands for the extrapolations referred to in the text. 
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excited mode of Hkink which for T/J < 1 closely follows the diluted kink picture already expected for equilibrium, 
i.e. pi — > 2/L. Also, preliminary evaluations of kink-kink correlations (ipi\n i rij\'i})i) show that these become negligible 
below T/J ~ 0.5, thus resembling the fully segregated equilibrium state. In passing it is worth pointing out that 
while Ai is doubly degenerate, the correlations and p\ remain indistinguishable in both of these eigenmodes. 
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FIG. 5: Density of kinks for the first excitation (solid line), and ground state (dashed line) of }lkink 
and J > 0. Lowermost dots indicate the equilibrium density (2/L) when T — > 0. 



Eq. JT8J ] with L — 24 



J < 



In studying numerically the AF dynamics, the approximants (|21[) pose severe difficulties at intermediate tempera- 
tures. Owing to low lying level crossings, now the spectrum gap is bisected in temperature regions having opposite 
monotonic behaviors. For instance, it is instructive to consider the case of L = 4 where this feature already appears 
in the spectrum of H S pin ■ By diagonalizing its associated 6x6 stochastic matrix, we find that there are two branches 
of eigenvalues crossing at K* = — j In 2, thus producing the non- analytic gap 



threefold degenerate if K > K* 
k , nondegenerate otherwise , 



(22) 



as seen in the uppermost curve of Fig. [6^. Alike this simple case, for larger sizes it is found numerically that Ai(L) 
remains degenerate only above certain temperatures below which level crossings occur, and Ai results in a nondegen- 
erate value. This is signaled by the emergence of the cusps observed in Fig.[6K. At LT regions the gap recovers the 
exponential decay referred to above but the data collapse now precludes finite-size estimations of effective exponents 
from Eq. (f2Tj) . Although level crossings come out at successively smaller temperatures, these can not keep pace with 
the increasing lattice sizes. To bypass these limitations we resorted to the antikink operator especially constructed 
for this region [Eq. (|20p ]. As can be seen in Fig.|BjD, now the LT spectrum no longer collapses so the approximants 
(|2"Tj) can be employed once more. However, as evidenced by the cusps of that figure, the problem of low lying level 
crossings yet persists. In nearing the crossing temperatures this brings about pronounced size effects in zl, partly 
because of the mismatches occurring between eigenvalue branches, alike those observed in Fig.|6K. 

All the above considerations result in fair convergent estimations of dynamic exponents both for LT and HT regimes, 
but as displayed in Fig. [71 the intermediate zone is out of reach. Despite this setback, what matters is the usefulness of 
the antikink representation (|20j) in the more interesting LT region which otherwise would have remained inaccessible. 
As in the F case, no extrapolations to the thermodynamic limit are needed here since already for T/\J\ < 0.3 our higher 
antikink approximants lie within the interval 1.98 ~ z%o ^ Z22 ^5 ^24 ~ 1.99, thus indicating a Glauber universality 
class for this ordering regime, cf. (loj . In the other extreme, when disorder prevails, already for T/\J\ > 3 the spin 
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FIG. 6: Spectrum gaps of (a) spin, and (b) antikink operators [Eqs. (110|) and (I20p respectively], for various inverse temperatures 
K — J/T, (J < 0). Sizes increase from top to bottom. In (a) the dashed line stands for the soluble case L = 4 already indicating 
the appearance of a cusp due to level crossings [see Eq. (l22|l ], a feature holding also for L = 6, 8, 24. At low temperatures 
Ai(L) decays as e _4 ' K ' with an almost size independent amplitude. In (b) this latter problem is remedied although level 
crossings for 8, 10, ... , and 24 sites now give rise to non-analiticities near which size effects become severe (cf. Fig. [7]). 
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FIG. 7: Effective dynamic exponents as defined in Eq. (f2Tjl . for Hamiltonians (fT0)l and ([20|) with J < and L — 10, 12, ... ,24 
increasing upwards. In nearing the cusps of Fig.[6]size effects become progressively pronounced, i. e. Ai (L) and Ai (L — 2) belong 
to different branches. High and low temperature regimes are both consistent with a typical diffusive exponent (224 — 1.99). 
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approximants derived from the data of Fig.[6£i converge within one percent to our LT exponents, thus suggesting for 
these latter a rather robust diffusive behavior. 

Finally, in Fig. [8] we exhibit the antikink density p\ for the first excited state of H a nti- As before, this confirms the 
diluted antikink picture expected at LT scales, in turn accounting for the convergent approximants obtained below 
T / J < 0.3. However, this time p\ does not quite follow there the vanishing equilibrium density of the plain AF vacuum 
state and adopts instead the limiting value referred to in Fig.[5j indicating the presence of two antikink excitations. 
At intermediate temperatures the jump of p\ just reflects the aforementioned level crossings. 

P 
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FIG. 8: Density of antikinks for the first level and ground state of operator (|20|l with L = 24 and J < (solid and dashed 
lines respectively). For comparison, horizontal dots indicate the equilibrium density of ferro kinks (J > 0) at T = 0, cf. Fig. [5] 
The discontinuity denoted by vertical dots stems from level crossings in the low lying spectrum of H a «ti (also, see cusp in the 
lowermost curve of Fig. [6b) . 



IV. CONCLUDING DISCUSSION 

To summarize, we have constructed symmetric representations of the Kawasaki dynamics in the Ising chain using 
a quantum spin analogy with both direct an dual processes. The kink description [7j], either in its F or AF versions 
[Eqs. (US]) and ]. is well suited to low temperature regimes where it is able to provide a rather fast finite-size 
convergence to dynamic exponents (Figs.HJa and [7]). For F couplings, the ordering kinetics arising from these latter 
(z ~ 3.11) tuns out to be slightly slower than the Lifshitz - Slyozov behavior [£(t) oc i 1 / 3 ] characteristic of higher 
dimensions [9(. By contrast, for AF interactions the kinetics is no longer activated by metastable states, and dynamic 
exponents converge rapidly (z ~ 1.99) to the typical diffusive values of the Glauber universality class [£(£) oc t 1 / 2 ]. 

On the other hand, the spin representation [ Eq. (|10[) ] is more appropriate for higher temperature scales where 
the kinetics remains slow and is still characterized by relaxation times diverging as L z (inset of Fig.HJi). For the 
F dynamics, the finite-size approximants of these exponents turn out to form a sequence of upper bounds which 
also extends down to low temperature regions (Fig. 0k), where it nicely complements the series of lower bounds 
emerging from the kink approach (Fig.lJb). On wider temperature scales, both representations are consistent with 
a nonuniversal set of dynamic exponents interpolating continuously between a subdiffusive ordering kinetics and the 
plain diffusive limit (Fig.[4j:). For J < 0, part of this intermediate region proved inaccessible due to size effects caused 
by low lying level crossings in the spectrum of both spin and antikink operators (cusps of Fig. [5]). However, as the 
dynamics of this case is not activated, the diffusive exponents obtained by both representations within their natural 
range of applicability (Fig. [7J suggest a rather universal (Glauber) dynamics throughout. 

The success of the kink description at low temperature regimes would probably allow to extend the ideas of this 
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work to other interesting Id-dynamics under instantaneous quenches (i.e. constant transition rates). For instance, 
the nonuniversal behavior alleged for the ordering kinetics of the alternating bond Kawasaki chain 1Q], might well be 
further investigated with our numerical approach. Other lines of research that would also be worth pursuing are certain 
dynamical processes involving composite objects (e.g . dimers) which exhibit strongly broken ergodicity as a result of 
having an extensive number of conservation laws [16( . As the issue of universality classes in nonequilibrium statistical 
systems is often linked to the existence of these latter [13], it would be important to determine whether dynamic 
exponents actually depend on the subspaces where the evolution takes place. Dual representations constructed 
appropriately for such processes could therefore provide reliable computations of these exponents within ordering 
regimes. Further work in this direction is in progress. 
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